EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-EP-2000-038 



10 March 2000 



Searches for Neutral Higgs Bosons 
in e"^e — Collisions 



Searches for neutral Higgs bosons in the Standard Model and the MSSM have 
been performed using data collected by the DELPHI experiment at a centre- 
of-mass energy of 188.7 GeV, corresponding to an integrated luminosity of 
158 pb _1 . These analyses are used, in combination with our results from 
lower energies, to set new 95% confidence level lower mass bounds on the Stan- 
dard Model Higgs boson (94.6 GeV/c 2 ) and on the lightest neutral scalar (82.6 
GeV/c 2 ) and neutral pseudoscalar (84.1 GeV/c 2 ) Higgs bosons in a represen- 
tative scan of the MSSM parameters. The results are also interpreted in the 
framework of a general two-Higgs doublet model. 



around 




DELPHI Collaboration 



Abstract 



(E. Phys. J. C17(2000)187/549) 



11 



P.AbreuH 
T.Allmendi: 
P.Andersso 
A.Augustm 
A Ra.rnn rnll ir 1 



A 



damE 



T.A. 



P. P. All 




P.AdzicH 
S.Almehedt 



A.AndreazzaO, S. AnHrinpaP 2 !. P.Antilag 
P.BambadeEi F.BaraoB, 




P.Baillo: 
M.Battagl, 
N.C.BenekQi 

M.BigiEl, M.ELBilenkyfcj, M-A.BizompM D.Blo 



K.Bel 



nil J 44 ! . 



M.Baubillierr 
A.C.Benvenutil 

n 



A.W.BomlandH, G.BQrisovE3, C^osicH, O.I 
iozzcfca, M. Brack 




IdiEa, N 
W-D.Ap. 
.Barbielli 
M.Begall 
M.Berggre 
H.M.BlonM 



T.Alderweirek 
apaneE 
Y.Arnoui 
R. Barbie] 
Behrmam 
D.Ber 
M.BonesinJ 



lekaei 



G.D.Ale 
AmatcO, 
B 

D.Y.Ba: 



snia 



R.Alemam 
E.G.Anassontz 
E.Augustin^ 
G.Barkeill! 



P. BeilliereS, Yu. Belokopytol 
D . BertrandEl M . Besancon _ 



M.Booneka: 



I.BozovicEj. 
B.Bnschhei 



uschmann 



E.Bo; 
P.Branchi 
S.Cabrer 




M.V. Castillo Gimenc: 



I.BoykcE 
T.Burarl 3 

F. Caren. 
Ph.Charpent^ 
V.Chorowi; 
D.Crennel 

G. Delia Riccal 4 ^, P.Delpierre^j, N.DemariaR A.De Angelis^, W.De Boertl C.De ClerccJ 2 ], B.De Lottd* 



,arrolH C.Ca: 
L.Cha uasa i 
ChudobaO, K.C 
, S.CrepeS, GX: 




tilk> 1 

ecchiaEJ , 
P.Collini 



G.A.Chel 




C.Bourdarioi 



arJ3, J-M.Brune 
T.CamporesjQ, 
V.ChabavLdl 
P.ChliapnikovEll 



J.Dolbe, 



GpEkspon, 




R.CpntriO, E.CortinaO, G.CopneES, F.Cossutt 
uevas Maestro^! S.Czellar fH J.DalrnauH, M. Davenport 



L.De Paul 
J-D.Duran< 
G.Fanoura] 

A. FirestoneEJ, LLElagmeye: 
F.Fulda-QuenzerEa, J.Fu: 
M.Gaspari2, U.Gasparin 

B. Golorjyila, G. Gom- 
el. Grahltl, E 
K.Hamachc: 
T.L.Hes; 
M.Hube: 
Ch. Jarlskq 




K.Dorob 
Eller 



M.DracosfcS 



A.Fer, 




J.DreesEl 



J-P.En_ 
-RibasEl, 



H.Foet i.Fokitisa, r E.Fontanellliil,^B.FranekE3, A. 1 



D.Ga: 



GazisEl, D.Gele&a, N. 
I.Gonzalez,_Caballer 
K.Grzel; 
V.HedberJJP 
S-O.Hol 
K.Hultqvis 
B.Jean-Marie ^ 



S.GamblinEa, M.G, 

hodbana?3 




elma 
I.Gi 
L.Go: 
S.Hah 
J. J.Hernande: 
S.Hoorelbekefel,, 
R. Jacobsso; 
E . K . Johansso 



G. Jarlski 
P. JuillotJlJ, Lujlingermannll 

B.P.KerRevara, YuTChokhloM B.A.KhomenkctZl, N.N.KhovanskM A.KiiskincrM_B.Kin; 



K.KarafasoulisB, S.KataanevasE3, E.C.Katsoufii 



H.Kleinfj, P.KLui 
Z.KrumsteinE-j 



0. KlappI 

E. Krizni 
J-P.Laugie 
J.Lem- 

1. LippJ 
L.Lyon: 

F. Mand: 
C.Martinez- 
C.Matteuzz 

G. Mc 
T.MoaNI, 

K.MuenictM M.MuldersI 
F.Naragh" 15 ! 
B.S.Nielsj 
A.Onofr 
J.Palacioi 
A.Passeril 

H. T.Philli 
P.Privite: 
A.L.Reai 
F.Richar< 
P.Rosinskyl 
A.Sadovaky 
F.ScuriEl, 
F.Simone 
T.Spass 
R.StrubE 
J. Thomas 
A.Tonazzi 



P.KokkiniasEl, V-KostioukhineM, C.Kou 



P.KubinedJ, 
:angast3, 
V.Lepeltie 
G.Loke: 



T.LesialM 
Lopesl 
J.R.MahoiiM, A.Mair^ 2 
R.Marcdl3, B.Man iiin 
nez-Vida]0. 




J.Kurowskal 
ederE3, EXedroitE. 



K.Kurvinen* 
cburc 3 



LcthuillicrEL J.LibbyEa, W.LiebigP 



J.M.Lop r - 

A. Male 




:oumelisa, 
I, J.W.Lami 
L.Leinone: 




El, P.S.L.Bootb!; 
~ >l) l T.J.y.Bowcock!; 

L.Bugge !f 
V.Canak- « 
M.Chapki-^ 
P.Chochulat 
H.B.Crawlevfe 
W.Da Silve^ 
|, A.De Miri' 
A.Duperrin ^ 
M.Eapirito Santo ^ 
S.Fichet 2 ± 
Fruhwirt v J 5 ^ 
C.Gasn 
R.Gokiclte 
V.GracccJii 
Hallgrer 
H.He: 
J.Hrubei 
I 1 ™ R.JanikE 



Dri 



F.Ferrdy, 



'rodesenQ, 
C.Gar 
F.Glegel 
il, Yu.Gou: 
S.Haidcrl 
P.Her 
M.Houlde 
P. Jaloch. 




F.MazziicatcEZI, 
WJ.Meyery : 
ellerEa, K.Moeni. 
.Mulet-Marqui 
F.L.Navarri. 
P.NiezirrawskE3, ^.Nikolenki 
R.Orava!^, G.OrazM K.Ostc, 
H.Palkaf4 Th.D.P,arjadopoulo 

IVLPegorar' 



s Marti i Garciafef^ J_Masi 



R.Lopez-FernE 
TX-M.Malmgrei 
J-C.Mari 
N.Mas 




D.Lo; 



MJMazzucatcEJ M.Mc Cubt 
iagko\C3, E/ 



M.R.Mong 
R.Muresan 




A.N. Si 



D.Sampsoni 
, A.IVL^ega: 




avasEil, K.Na 
V.Nomoko: 
A.Ourao 
K.Papag> 
L.Peralt. 



lioreH, L 
X.Moreaul 
W.J.Mur 
,ock 



B.Mupipta, G.Myat 
N.Neuf. 
V.Obrazlso 




T.PodobnikH, M 

n 



Ragazz 
hnt!l, D.R 
O.Rohnd^fl 



P. Jonsson ^ , C . Jorar 
Keranenl^, G. Kernel 
A.Kip^igN, N.J.Kjae 
Krammei 52 
V.Lapin^ 
R.Leitner ^ 
AJLipniackau ti 
,sE3, P.Lut2 
V.Malychev 

A.Markou l 2 
F.Matorras * 
R.Mc Nultyp. 
joernmar 
U.Muellei£i 
T.Myklebusipi 
R.Nicolaidor 
Olshevsk 
R.PainE3, R.PaivaE; 
F.Parodjii U.Parzefalf; 
C.Petridm£3, A.Petrolin 
PoropatC _^V.Pozdniakov l 2. 
J.Ramesy, P.N.Ratoffc 
L.K.Resvanis 



3, S.MaltezosQ, 
3, C.MariotJ 
iyiannopoulos 
R.Mc K| 
W.A_Mitarofl 5: 




k 2 l 



P. B. Rent, 
RoncheS' 



I.RosenbergE 



.Ruhlmarin-KlcidciEl, AjRuizC 2 ], H.Saarikkot3, Y.Sacquini^ 
SanniflcO, Ph . Srliwcm 1 i n pp 4 l . B Sr.hwerin pJ 54 ! . U^chwickerath _ 
Seiherty, R.Sekuliry, R.C.Sheilardl, 



M.Sieh 



E.SpiritiEJ, S . Squarci al 14 l Q . Stan esr.i I 40 !. S Sta.ni- 



.R.Smithri. QjSoIovianovCj, A . S opcza.l 



L.Simarc i 1 
„ Sosnowsk*" 



M.Sta 



tzkita, K.StevensonE3, A.StocchO, J.Strausf^ 



B.StuguH, M.SzczekowskiH, M.SzeptyckaEJ, JJ.TabarellEj A.TaffarcEa, F.TegenfeldtH, E.Terranova^ 
til, L.G.TkatclyevtZI, M.TohinB, S.Todorovay, A.TomaradzeB, B.Tomc^ 



J. Tim merman; 

L.TortoraH P.TortosaH, G.Transtromeit3, D.TreilleH, G.TristramH, M.TrochimczukH, C.TronconEf 



I.Tmi 
sa0, 




I.A.Tyapk 
W.Van 
G.Veg: 
Vitalc 



riasEl, 



M-L.Turlue 
P. Van Daml 
I.Van Vulp 
D.Vilanov; 
A.J.Washbroo. 
O.Yushchenki 

Ph.Zolleita, G.C.ZucchelliHl G.Zumerld2Zl 



b. . .zamarias H, O.Ullalani 



BoeckEl, 
L.Ventu: 
E.Vlaso 




Van Doni) 
W. Venus 



V.Uvarov0 n i 



.ValentI 



E.Vallazz. 



TP 



C.Vander Velde: , 



J. Van EldikH'c3, A. Van LysebettenEl, N.van Remortelt , 
RVerbeurekl, MA^erlatcO, _L.S .VertogradcpM V.Verz: 3 1 , 



A.S.VodopyanovtZl, G.VoulgarisH, VVrbaM, H.WahknEj, CWalck 



G.WolfQ, J.^ 
A.Zintchonkol 17 ' 



C.Weisei0,_D. Wickets, J.ILWickensB, G.ILWilkinsox£!l, M.WinterE3, M.Wij 
A.Zalewskay, P.ZaWskiE], D.ZavrtaniklfJ, E.ZevgolatakosEI , N.I.Zimint 



1 Department of Physics and Astronomy, Iowa State University, Ames IA 50011-3160, USA 
2 Physics Department, Univ. Instelling Antwerpen, Universiteitsplein 1, B-2610 Antwerpen, Belgium 
and IIHE, ULB-VUB, Pleinlaan 2, B-1050 Brussels, Belgium 

and Faculte des Sciences, Univ. de l'Etat Mons, Av. Maistriau 19, B-7000 Mons, Belgium 
3 Physics Laboratory, University of Athens, Solonos Str. 104, GR-10680 Athens, Greece 
4 Department of Physics, University of Bergen, Allegaten 55, NO-5007 Bergen, Norway 
5 Dipartimento di Fisica, Universita di Bologna and INFN, Via Irnerio 46, IT-40126 Bologna, Italy 
6 Centro Brasileiro de Pcsquisas Fi'sicas, rua Xavier Sigaud 150, BR-22290 Rio de Janeiro, Brazil 
and Depto. de Fisica, Pont. Univ. Catolica, CP. 38071 BR-22453 Rio de Janeiro, Brazil 
and Inst, de Fisica, Univ. Estadual do Rio de Janeiro, rua Sao Francisco Xavier 524, Rio de Janeiro, Brazil 
7 Comenius University, Faculty of Mathematics and Physics, Mlynska Dolina, SK-84215 Bratislava, Slovakia 
8 College de France, Lab. de Physique Corpusculaire, IN2P3-CNRS, FR-75231 Paris Cedex 05, France 
9 CERN, CH-1211 Geneva 23, Switzerland 

10 Institut de Recherches Subatomiques, IN2P3 - CNRS/ULP - BP20, FR-67037 Strasbourg Cedex, France 
11 Now at DESY-Zeuthen, Platanenallee 6, D-15735 Zeuthen, Germany 

12 Institute of Nuclear Physics, N.C.S.R. Demokritos, P.O. Box 60228, GR-15310 Athens, Greece 

13 FZU, Inst, of Phys. of the C.A.S. High Energy Physics Division, Na Slovance 2, CZ-180 40, Praha 8, Czech Republic 
14 Dipartimento di Fisica, Universita di Genova and INFN, Via Dodecaneso 33, IT-16146 Genova, Italy 
15 Institut des Sciences Nucleaires, IN2P3-CNRS, University de Grenoble 1, FR-38026 Grenoble Cedex, France 
"Helsinki Institute of Physics, HIP, P.O. Box 9, FI-00014 Helsinki, Finland 

17 Joint Institute for Nuclear Research, Dubna, Head Post Office, P.O. Box 79, RU-101 000 Moscow, Russian Federation 

18 Institut fur Experimentelle Kernphysik, Universitat Karlsruhe, Postfach 6980, DE-76128 Karlsruhe, Germany 

19 Institute of Nuclear Physics and University of Mining and Metalurgy, Ul. Kawiory 26a, PL-30055 Krakow, Poland 

20 Universite de Paris-Sud, Lab. de l'Accelerateur Lineaire, IN2P3-CNRS, Bat. 200, FR-91405 Orsay Cedex, France 

21 School of Physics and Chemistry, University of Lancaster, Lancaster LAI 4YB, UK 

22 LIP, 1ST, FCUL - Av. Elias Garcia, 14-1°, PT-1000 Lisboa Codex, Portugal 

23 Department of Physics, University of Liverpool, P.O. Box 147, Liverpool L69 3BX, UK 

24 LPNHE, IN2P3-CNRS, Univ. Paris VI et VII, Tour 33 (RdC), 4 place Jussieu, FR-75252 Paris Cedex 05, France 

25 Department of Physics, University of Lund, Solvcgatan 14, SE-223 63 Lund, Sweden 

26 Universite Claude Bernard de Lyon, IPNL, IN2P3-CNRS, FR-69622 Villeurbanne Cedex, France 

27 Univ. d'Aix - Marseille II - CPP, IN2P3-CNRS, FR-13288 Marseille Cedex 09, France 

28 Dipartimento di Fisica, Universita di Milano and INFN, Via Celoria 16, IT-20133 Milan, Italy 

29 Universita degli Studi di Milano - Bicocca, Via Emanuelli 15, IT-20126 Milan, Italy 

30 Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

31 IPNP of MFF, Charles Univ., Areal MFF, V Holesovickach 2, CZ-180 00, Praha 8, Czech Republic 
32 NIKHEF, Postbus 41882, NL-1009 DB Amsterdam, The Netherlands 

33 National Technical University, Physics Department, Zografou Campus, GR-15773 Athens, Greece 
34 Physics Department, University of Oslo, Blindern, NO-1000 Oslo 3, Norway 
35 Dpto. Fisica, Univ. Oviedo, Avda. Calvo Sotelo s/n, ES-33007 Oviedo, Spain 
36 Department of Physics, University of Oxford, Keble Road, Oxford OX1 3RH, UK 
37 Dipartimento di Fisica, Universita di Padova and INFN, Via Marzolo 8, IT-35131 Padua, Italy 
38 Rutherford Appleton Laboratory, Chilton, Didcot OX11 OQX, UK 

39 Dipartimento di Fisica, Universita di Roma II and INFN, Tor Vergata, IT-00173 Rome, Italy 

40 Dipartimento di Fisica, Universita di Roma III and INFN, Via della Vasca Navale 84, IT-00146 Rome, Italy 

41 DAPNIA/Service de Physique des Particules, CEA-Saclay, FR-91191 Gif-sur-Yvette Cedex, France 

42 Instituto de Fisica de Cantabria (CSIC-UC), Avda. los Castros s/n, ES-39006 Santander, Spain 

43 Dipartimento di Fisica, Universita degli Studi di Roma La Sapienza, Piazzale Aldo Moro 2, IT-00185 Rome, Italy 

44 Inst, for High Energy Physics, Serpukov P.O. Box 35, Protvino, (Moscow Region), Russian Federation 

45 J. Stefan Institute, Jamova 39, SI-1000 Ljubljana, Slovenia and Laboratory for Astroparticle Physics, 

Nova Gorica Polytechnic, Kostanjcviska 16a, SI-5000 Nova Gorica, Slovenia, 

and Department of Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia 
46 Fysikum, Stockholm University, Box 6730, SE-113 85 Stockholm, Sweden 

47 Dipartimento di Fisica Spcrimentale, Universita di Torino and INFN, Via P. Giuria 1, IT-10125 Turin, Italy 
48 Dipartimento di Fisica, Universita di Trieste and INFN, Via A. Valerio 2, IT-34127 Trieste, Italy 

and Istituto di Fisica, Universita di Udine, IT-33100 Udine, Italy 
49 Univ. Federal do Rio de Janeiro, CP. 68528 Cidade Univ., Ilha do Fundao BR-21945-970 Rio de Janeiro, Brazil 
50 Department of Radiation Sciences, University of Uppsala, P.O. Box 535, SE-751 21 Uppsala, Sweden 
51 IFIC, Valencia-CSIC, and D.F.A.M.N., U. de Valencia, Avda. Dr. Moliner 50, ES-46100 Burjassot (Valencia), Spain 
52 Institut fur Hochenergiephysik, Osterr. Akad. d. Wissensch., Nikolsdorfergasse 18, AT-1050 Vienna, Austria 
53 Inst. Nuclear Studies and University of Warsaw, Ul. Hoza 69, PL-00681 Warsaw, Poland 
54 Fachbereich Physik, University of Wuppertal, Postfach 100 127, DE-42097 Wuppertal, Germany 



1 



1 Introduction 

In the framework of the Standard Model (SM) there is one physical Higgs boson, 
H, which is a neutral CP-even scalar. At LEP II the most likely production process is 
through the s-channel, e + e~— > Z* ^HZ. The W + W~ and ZZ fusion t-channel production 
processes in some of the channels described below are not considered here, but their 
contribution is typically below 10% in the range of masses considered in this study. 

The results of the search for the SM Higgs are also interpreted in terms of the lightest 
scalar Higgs boson, h, in the Minimal Super-symmetric Standard Model (MSSM). This 
model predicts also a CP-odd pseudo-scalar, A, produced mostly in the e + e~— > hA process 
at LEP II. This associated production is also considered in this paper. 

With the data taken previously at y/s= 183 GeV DELPHI excluded a SM Higgs 
boson with mass less than 85.7 GeV/c 2 |l|], and set limits on h and A of the MSSM of 
74.4 GeV/c 2 and 75.3 GeV/c 2 respectively. The present analyses therefore concentrate on 
masses between these and the kinematic limit. Note that the LEP Higgs working group 
has found mass limits of 89.7 GeV/c 2 for H, 80.1 GeV/c 2 for h and 80.6 GeV/c 2 for A, 
under assumptions generally referred to as the benchmark scan, when combining the data 
of the four LEP experiments from data taken up to 183 GeV. 

In the HZ channel, all known decays of the Z boson have been taken into account 
(hadrons, charged leptons and neutrinos) while the analyses have been optimised either 
for decays of the Higgs into bb, making use of the expected high branching fraction 
of this mode, or for Higgs boson decays into a pair of r's. A dedicated search for the 
invisible Higgs boson decay modes will be reported separately. The hA production has 
been searched for in the 4b and bb r + r~ channels. 

There are separate analyses for the different decay modes of the Higgs and Z bosons. 
Some common features are discussed in Sect. |], the H/i + yU _ and He + e~ channels in Sect. £|, 
channels involving jets and r's in Sect. |5] and ~Rvv in Sect. || Purely hadronic final states 
are discussed in Sect. [f[ The results are presented in Sect. §. 

2 Data samples overview and the DELPHI detector 

For most of the data collected in 1998, LEP was running at energies around 189 GeV. 
DELPHI recorded an integrated luminosity of (158±1) pb _1 at a mean energy of 
188.7 GeV. 

Large numbers of background and signal events have been produced by Monte Carlo 
simulation using the DELPHI detector simulation program ||. The size of these sam- 
ples is typically about 100 times the luminosity of the collected data. Background was 
generated with PYTHIA § and K0RALZ [|] for (e+e - -»• ff 7), PYTHIA and EXCALIBUR 
for the four-fermion background and TW0GAM |7j and BDK || for two-photon processes. 
BABAMC was used to simulate Bhabha events in the main acceptance region. 

Signal events were produced using the HZHA [JTOfl generator. For the SM process the 
Higgs mass was varied in 5 GeV/c 2 steps from 70 GeV/c 2 to 100 GeV/c 2 , while for hA , 
the A mass was varied between 70 and 90 GeV/c 2 with tan/? (the ratio of the vacuum 
expectation values of the two doublets) either 2 or 20. This fixes the h mass, almost 
equal to for tan/3 = 20 and significantly lower than m\ if tan/3 = 2. 

The HZ simulated samples were classified according to the Higgs and Z boson decay 
modes. For He + e~, H/i + /i~ and Hi/v the natural SM mix of H decay modes was permitted. 
In the Hqq channel the rr decay mode was removed, and we generated separately the 
two channels involving r leptons for which one of the bosons is forced to decay to r's and 
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the other hadronically. Finally, for the hA simulations final states involving either four 
b quarks or two b quarks and two r's were simulated. Efficiencies were defined relative 
to these states. The size of these samples varied from 2000 to 20,000 events. 

The detector was unchanged from the previous data taking period. Thus we refer 
to our previous publication JIJ for a short description. More details can be found in 
references P,|TT] . 



3 Common features for all channels 

3.1 Particle selection 

In all analyses, charged particles are selected if their momentum is greater than 
100 MeV/c and if they originate from the interaction region (within 10 cm along the 
beam direction and within 4 cm in the transverse plane). Neutral particles are defined 
either as energy clusters in the calorimeters not associated to charged particle tracks, or 
as reconstructed vertices of photon conversions, interactions of neutral hadrons or de- 
cays of neutral particles in the tracking volume. All neutral clusters of energy greater 
than 200 MeV (electromagnetic) or 500 MeV (hadronic) are used; clusters in the range 
100-500 MeV are considered with specific quality criteria in some analyses. The ir^ mass 
is used for all charged particles except identified leptons, while zero mass is used for 
electromagnetic clusters and the K° mass is assigned to neutral hadronic clusters. 



3.2 b-quark identification 



The method of separation of b quarks from other flavours is described in ||12|| , where 
the various differences between B-hadrons and other particles are accumulated in a sin- 
gle variable, hereafter denoted x\> for an event and x l h for jet i. A major input to the 
combined variable is the probability that all tracks in a group originate from the inter- 
action point. Xb combines this probability with information from secondary vertices (the 
mass computed from the particles assigned to the secondary vertex, the rapidity of those 
particles, and the fraction of the jet momentum carried by them) and also the trans- 
verse momentum (with respect to its jet axis) of the leptons, using the likelihood ratio 
technique. Increasing values of x^ correspond to increasingly 'b-like' events (or jets). 

The procedure is calibrated on events recorded in the same experimental conditions 



at the Z resonance. The performance of the combined b-tagging is described in Ref. |13 



and the impact parameter tagging in Ref. ||14|| . The overall performance is illustrated in 
Figure [p. 

A careful study of possible systematic effects, including data versus simulation agree- 
ment at the Z pole (checked inclusively, per flavour and for multi-jet events), leads to an 
overall relative b-tagging uncertainty below 5%, varying slightly with the exact tagging 
value used. At high energy, an inclusive comparison of data with simulation confirms this 
number. The gluon splitting rates into bb and cc have been rescaled in the simulation 
according to the DELPHI |15| measurement. In addition, a 50% uncertainty on these 



splitting rates is applied to the qq (7) background estimate. 
3.3 Constrained fits 



In most channels a constrained fit [0 is performed to extract the Higgs mass, and often 
to reject background processes as well. If only total energy and momentum conservation 
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are imposed then the fit is referred to as '4-C, while some fits require the Z mass as well, 
either as a fixed value, or taking into account the Breit-Wigner shape of the Z resonance. 
In both cases such fits are referred to as '5-C In order to allow the removal of most 
of the radiative return to the Zevents, an algorithm has been developed []T7| in order to 
estimate the effective energy of the e + e~ collision. This algorithm makes use of a '3-C 
kinematic fit in order to test the presence of an initial state photon parallel to the beam 
direction and then lost in the beam pipe. This effective centre-of-mass energy is called 
\fs i throughout the paper. 



3.4 Confidence level calculations 

The procedure used to compute the confidence levels is the same as that used in our 
previous publications [0,0 but the discriminant information is now two-dimensional: 
the first variable used is the reconstructed mass, the second one is either the x\, (for the 
electron and muon channels) or the likelihood (for all other channels). As far as the mass 
information is concerned, the reconstructed Higgs boson mass is used in the hZ channels 
and the sum of the reconstructed h and A masses in the hA channels (for the pairing 
with minimal mass difference in the four-jet channel). In order to make full use of the 
information contained in the second variable the selections are looser than in the past: 
the method used for deriving the confidence levels ensures that adding regions of lower 
signal and higher background can only enhance the performance expected from a tighter 
selection. Since the distributions, represented as two-dimensional histograms, are derived 
from simulation samples, the limited statistics in some bins are a potential problem: 
statistical fluctuations can artificially increase the expected sensitivity. This possible 
systematic shift of the confidence level has been estimated comparing the expected results 
using the full simulation sample with those derived from fractions of this sample. The 
bin sizes were carefully chosen to keep full sensitivity while avoiding any significant bias 
caused by this effect. 

4 Higgs boson searches in events with jets and elec- 
trons or muons 



4.1 Electron channel 

The analysis follows what has already been published |l| , with the following improve- 
ments in the selection cuts. To reinforce the Bhabha veto, the preselection described 
in [0] has been complemented by a rejection of electron candidate pairs having acopla- 
narity (defined as the supplement of the angle between the transverse momenta of the 
two electrons) below 3 degrees and energies above 40 GeV. To allow for the tau decays 
of the Higgs boson while keeping a good purity, the requirement on the minimum event 
charged multiplicity has been raised to 8 except if the recoiling system from the electron 
candidate pair is made of two jets, each with a charged multiplicity lower than or equal 
to 3 and with a mass below 2 GeV/c 2 . This defines the preselection. 

The energy of the slower/faster electron is required to be above 15/20 GeV. Electron 
isolation angles with respect to the closest jet are required to be more than 20° for 
the most isolated electron and more than 8° for the other. Global kinematic fits |TE[] 
are performed, imposing total energy and momentum conservation and constraining the 
invariant mass of the e + e~ system to mi (5-C fits). If the fit probability is below 10 -8 , the 
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fit procedure is redone for fixed values of the e + e~ mass between 36.5 and 105 GeV/c 2 , 
in order to allow for the tails of the Z mass distribution. For each mass a combined 
probability is determined as the product of the \ 2 probability times the Breit-Wigner 
probability for the Z mass, and the mass giving the maximum combined probability is 
retained. Events with a combined probability below 10~ 8 are rejected. As the search is 
restricted to high mass Higgs bosons produced in association with a Z particle, the sum 
of the masses of the electron pair and of the recoiling system as given by the kinematic fit 
is required to be above 150 GeV/c 2 and the difference in the range from -100 GeV/c 2 to 
50 GeV/c 2 . The fitted recoil mass and the global b-tagging variable Xb are used for the 
two-dimensional calculation of the confidence levels. 

Table [l] shows the effect of the cuts on data, simulated background and signal events. 
The agreement between data and background simulation is illustrated at preselection 
level in Figure [| which shows the distributions of the main analysis variables, namely, 
the slow electron energy, the fitted mass of the jet system, the minimum electron isolation 
angle and the global event b-tagging variable x^. 

The final background amounts to 6.63 ±0.26 (stat-)-a93 (syst.) events, and is mainly 
due to e + e~qq (ZZ) events. Illustrations of the two-dimensional distribution, used as 
input for the confidence level computations are shown in figure |3] for data, simulated 
background and signal events (for mn =95 GeV/c 2 ). Table shows the selection efficiency 
for different Higgs boson masses. The systematic uncertainties have been evaluated as 
described in [ffj. Among the events selected in data, one has a high Xb value and thus is 
kept in the final mass-plot, for which the supplementary cut Xb > —1.8 was applied. 

4.2 Muon channel 

The analysis is based upon the same discriminant variables as in [p]], but the selection 
criteria have been re-optimised, as explained in 0, to afford the best sensitivity to the 
expected signal at 188.7 GeV. The preselection remains unchanged except that events 
must now have at least nine charged particle tracks. Two muons are required with 
opposite charges, momenta greater than 34 GeV/c and 19 GeV/c, and an opening angle 
larger than 10°. The muon isolation angle with respect to the closer jet must be greater 
than 16° for the most isolated muon and greater than 8° for the other one. A 5-C 
kinematic fit is then performed to test the compatibility of the di-muon mass with the Z 
mass. Events are kept only if the fit converges. The fitted recoiling mass is chosen as the 
first discriminant variable for the two-dimensional calculation of the confidence levels. 
The second variable is the global b-quark variable X\>. 

Table [I] details the effect of the selections on data and simulated samples of background 
and signal events. The agreement of simulation with data is good. This can also be seen 
in Figure |], which shows the total energy of the charged particles, the momentum of the 
fast muon candidate and the content in b-quark of the event after the preselection. The 
isolation of the muons with respect to the closest jet is also given after the lepton pair 
selection. At the end of the analysis, 5 events are selected in the data in good agreement 
with the expected background of 5.09 ± 0.19(stat.) ±0.21(syst.) events coming mainly 
from ZZ. Finally the signal efficiencies for different Higgs boson masses are given in 
Table |^. The systematic uncertainties on background and efficiencies have been derived 
as explained in The two events with the largest values of x^ are kept for the final 
mass-plot. 
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5 Higgs boson searches in events with jets and taus 

Three channels are covered by this analysis, two for the SM, depending of which boson 
decays into r + r~, and one for the MSSM. Hadronic events are selected by requiring at 
least ten charged particles, a total reconstructed energy greater than OAy/s , a recon- 
structed charged energy above 0.2y/s and an effective centre-of-mass energy, a/s 7 , greater 
than 120 GeV. This defines the preselection. 

A search for r lepton candidates is then performed using a likelihood technique. Clus- 
ters of one or three charged particles are first preselected if they are isolated from all other 
particles by more than 10°, if the cluster momentum is above 2 GeV/c and if all parti- 
cles in a 10° cone around the cluster direction make an invariant mass below 2 GeV/c 2 . 
The likelihood variable is calculated for the preselected clusters using distributions of the 
cluster momentum, of its isolation angle and of the probability that the tracks forming 
the cluster come from the primary vertex. Pairs of r candidates with opposite charges 
and an opening angle of at least 90° are selected using a cut on the product of their 
likelihoods, considering both the 1-1 (where the two r leptons decay to one prong) and 
1-3 (with at least one r decaying to three prongs) topologies. As an example, Figure [|a 
shows the r selection likelihood distribution for the selected events in the 1-1 topology. 

Two slim jets are then reconstructed with all particles (charged and neutral) inside a 
10° cone around the cluster directions. The rest of the event is forced into two jets using 
the DURHAM algorithm. The slim jets are constrained to be in the 20°< 9 T < 160° 
polar angle region to reduce the Ze + e~ background, while the hadronic dijet invariant 
mass is required to be between 20 and 110 GeV/c 2 in order to reduce the qq(7) and 
Z7* backgrounds. The jet energies and masses are then rescaled, imposing energy and 
momentum conservation, in order to improve the estimation of the masses of the dijets 
(t + t~ and qq ). Both dijets are required to have a rescaled mass above 20 GeV/c 2 and 
below y/s , and each hadronic jet must have a rescaling factor in the range 0.4 to 1.5. 

The remaining background comes from genuine £ + £~qq events. In order to reject the 
qq e + e~ and qq yU + /i~ backgrounds the measured mass of the leptonic dijet is required to 
be between 10 and 80 GeV/c 2 and its electromagnetic energy to be below 60 GeV (see 
Figure |5]c). The effect of the selections on data, simulated background and signal events 
is given in Table [3], while the selection efficiencies are summarised in Tables || (SM) 
and |3] (hA). Systematic uncertainties have been estimated by moving each selection cut 
according to the resolution of the corresponding variable. The main contributions are due 
to the r + T _ invariant mass and electromagnetic energy. The total systematic uncertainties 
amount to ±6% on signal efficiencies and ±11% on the background. 

At the end of the above selections, the reconstructed Higgs boson mass is estimated 
from the sum of the rescaled dijet masses in the hA channel and by subtracting the 
nominal Z mass in the hZ channels (Figure ^e). Besides this reconstructed mass, the two- 
dimensional calculation of the confidence levels makes use of a discriminating variable, 
again using a likelihood technique. This variable is built from the distributions of the 
rescaling factors of the r jets, the r momenta and the global b-tagging x\> variable (see 
Figure |5|g). Since the three possible r + r~qq signals are analysed in the same way, the 
confidence level computation uses only one global r + r~ qq channel. At each test point, 
the signal expectation and distribution in this channel are obtained by summing the 
contributions from the three signals weighted by their expected rates. 
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6 Higgs boson searches in events with jets and miss- 
ing energy 

The analysis starts with a preselection which is done in two steps. The first step aims 
at reducing the 77 contamination and requires a total charged multiplicity greater than 
10 (with at least one charged particle with a transverse momentum above 2 GeV/c), a 
total charged energy greater than 30 GeV and the sum of the transverse energies of the 
charged particles with respect to the beam axis greater than 28 GeV. The total transverse 
momentum has to be greater than 2 GeV/c. Furthermore, events where both the total 
transverse momentum and the largest single charged particle transverse momentum are 
less than 5 GeV/c have also been rejected. After these cuts, the 77 contamination is 
reduced to 1.5% of the total background, which is now dominated by 99(7) events. 



Then, jets are reconstructed from the event particles using the LUCLUS |L9| algorithm 
with the DURHAM distance (y cu t = 0.005). The results will be hereafter referred to as 
"free-jet clustering" . Events are also forced in a two-jet topology using the same algorithm 
(with a result referred to as "two-jet clustering") in order to tag specifically qq(j) events 
with the photon emitted along the beam axis. The rest of the preselection is designed to 
remove a large fraction of the remaining background without affecting the signal efficiency 
using the good discrimination between background and signal in the distributions of 
many analysis variables. The selection requires the most energetic electromagnetic cluster 
associated to a charged particle to be lower than 25 GeV (or 10 GeV if the charged particle 
associated to the cluster failed the charged particle selection criteria), the effective centre- 
of-mass energy, a/s 7 , to be greater than 100 GeV, the absolute value of the cosine of the 
polar angle of the missing momentum to be lower than 0.98, the most forward jet in the 
free-jet clustering to be more than 16° from the beam axis, the fraction of electromagnetic 
energy per jet in the free-jet clustering to be lower than 0.8, the energy deposited in the 
forward region within 30° around the beam axis to be smaller than 20 GeV and the energy 
of the more (less) energetic jet in the two-jet clustering to be between 30 and 90 GeV (10 
and 60 GeV). This defines the preselection. 

The final discrimination between signal and background is achieved through a multi- 
dimensional variable built using the likelihood ratio method. The input variables are the 
effective centre-of-mass energy y/s', the global b-tagging variable x^, the missing momen- 
tum P m is, and the cosine of its polar angle, the charged multiplicities of the jets in the 
free-jet clustering, the energy of the most energetic jet in the two-jet clustering, Ej et i, 
the acoplanarity, defined as the supplement of the angle between the transverse momenta 
of the two jets in the two-jet clustering, the maximal (over all particles in an event) 
transverse momentum with respect to the axis of the closest jet in the two-jet clustering, 
Ptmaxi and the output of a veto algorithm based on the response of the lead-scintillator 
counters installed at polar angles around 40° to detect photons crossing this insensitive 
region of the electromagnetic calorimeters. 

The distributions at preselection level of some of the input variables are shown 
in Figure |B| while that of the discriminant variable is given in Figure [?J Figure |^ 
shows the expected background rate as a function of the efficiency for a Higgs mass 
mh=95 GeV/c 2 when varying the cut on the discriminant variable. The final selection 
yields a total background of 27.8 ± 1.0(stat.) at 54.3% efficiency for the signal. The 
number of observed events is 27. Table |l] details the effects of the selections on data and 
simulated samples of background and signal events while the selection efficiencies as a 
function of the Higgs boson mass are summarised in Table |[ 
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The systematic uncertainty is dominated by the imperfect modelling of the energy 
flow. The corresponding error has been estimated by comparing data and simulation in 
test samples of Z7 events at high energy and data taken at the peak. It amounts to 11.0% 
relative to the background. Other sources include the imperfect modelling of b-tagging 
and jet angular resolutions, the dependence on the jet algorithm and uncertainties in 
cross-sections. This yields a total systematic error of ±12.4%. 



7 Higgs boson searches in pure hadronic events 

The aim of the four-jet preselection is to eliminate radiative and 77 events and to 
reduce the QCD and Z7* background. This preselection, common to HZ and hA analyses, 
has not changed with respect to last year's analysis , except that the number of charged 
particles for the di-jet recognised as the Z is required to be greater than or equal to five, 
in order to remove events in which the Z decays into charged leptons. 



7.1 The HZ four-jet channel 

The present analysis is an update of the method used by DELPHI at 183 GeV |T|. 
Events are selected using a discriminant variable which is defined as the ratio of like- 
lihood products for signal and background hypotheses for a set of quantities having a 
different behaviour in the two cases. These variables can be divided into two categories 
related respectively to the shape and to the b-content of the events. The six shape vari- 
ables are those defined previously and a new quantity: the fitted mass of the di-jet 
assigned to the Z. The agreement between data and background simulation is illustrated 
at preselection level in Figure |9| which shows the distributions of four analysis variables. 

In the previous version of this analysis the Z boson mass was fixed at its nominal 
value and the "best" pairing of jets to select the Higgs and Z candidates was found 
by maximising an expression in which the b-content of the different jets and the \ 2 
probability of the five-constraint fit contribute. When the production of the Higgs boson 
is close to the kinematic limit, the Z mass distribution no longer has a Breit-Wigner 
shape centred on mz ■ The previous procedure has been generalised using the Z mass 
distribution given by the simulation for a fixed Higgs mass equal to 95 GeV/c 2 . For values 
of the Higgs mass which differ from 95 GeV/c 2 , the same Z mass distribution has been 
used, even if not optimal, in order to be independent of the assumed value for the Higgs 
mass. The resulting mass distribution, for preselected events in data and simulation, is 
shown in Figure ITD. 



Events originating from the signal and from the background are separated using the 
value of a discriminant variable. This variable combines the information from shape 
variables to reduce the QCD background and from the b-tagging variable to reduce the 
contribution of W pairs. The most effective variable against the W + W~ background 
was found to be x l b , the combined b-tagging variable |12| measured for each jet. The 



likelihoods that each event is of HZ, W + W or QCD origin are evaluated. The final 
discriminant variable is obtained as the ratio between signal and background likelihoods. 



Figure |TT| shows the number of expected Monte-Carlo and observed data events, as a 
function of the efficiency on the signal when cutting on this discriminant variable. 

The two-dimensional distribution obtained by combining the Higgs mass estimate and 
the likelihood ratio is used for the final limit calculation. The most background-like 
events are suppressed by demanding that the (log) likelihood ratio be greater than -1.0. 
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The systematic uncertainties have been evaluated by considering a b-tagging, a QCD 
related and a 4-fermion part, independently, together with an uncertainty on the cross- 
sections of all processes, resulting in an overall relative uncertainty of 7.5% at the final 
selection level. 



7.2 The hA four— b channel 

A likelihood method has been applied to search for hA production in the four-jet chan- 
nel. After the common four-jet preselection, tighter cuts were applied to the remaining 
events, namely, a cut in the parameter of the DURHAM algorithm (y cu t > 0.003) is 
imposed as well as the requirement of at least two charged particles per jet. Finally, an 
event is rejected if its maximum inter-jet energy difference is greater than 70 GeV. The 
resulting number of expected events and the signal efficiencies after this preselection are 
given in Table |l|. 

The following eight variables are combined in the likelihood: the event thrust, the 
second and fourth Fox- Wolfram moments, H2 and H±, the minimal (among the three 
possible pairings of jets) di-jet-masses difference, the production angle of the candidate 
bosons, the sum (over the four jets) of the b-tag jet variable, the minimum di-jet b- 
tag and the number of secondary vertices. For each event, the measured value of each 
of these discriminant variables is compared with probability density functions obtained 
from simulated events. 

Figure [12] shows the resulting efficiency versus the total background varying the like- 
lihood cut (assuming mj^ = 80 GeV/c 2 and tan/3 = 20). The final cut value is chosen 
depending on the efficiency-background point desired; for the derivation of the limits, a 
cut on the likelihood output at 2.0 has been chosen. A total of 13 events is observed while 
11.41 are expected. A more stringent cut on the likelihood output corresponding to a 



requirement at 3.05 is used for the mass plot (see Fig [13]). This yields a total background 
of 4.72 ± 0.22 coming from qq(7)( 2.11 ± 0.14) and 4-fermion processes (2.61 ± 0.17) and 
3 events selected in data with a sum of their di-jet masses of 147, 180 and 176 GeV/c 2 . 
Efficiencies obtained for different masses and tan (3 are summarised in Table 0. 

To check systematic uncertainties on the total background due to the modeling of 
the shape of the probability density functions, the training and validation sample were 
exchanged and the analysis repeated. The uncertainty on the total background due to 
this effect has been estimated at the level of 5% and has been added (quadratically) to 
the other sources of errors, in particular the one coming from the b-tagging estimation, 
resulting in a final relative systematic uncertainty of 8%. 



8 Summary and results 

The results of the searches presented in the previous sections can be translated into 
exclusion limits on the masses of the neutral Higgs bosons in the SM and MSSM. 

8.1 Summary 

For each analysis of the HZ and hA channels at 188.7 GeV, the integrated luminosity, 
the expected backgrounds and their errors, and the number of observed events at various 
levels of the analyses are summarised in Table [3]. Within each channel, the penultimate 
line represents the inputs for the confidence level calculations ("final selection"), while 
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Selection 


Data lotal Q.Q.VTJ 


4 fermion Efficiency 




background 








Electron channel 155.4 pb 


— i 




Preselection 


1290 1227.4 924 


267 


78.5 


tight lepton id. 


28 27.9± 0.9 13 


12.4 


61.0 


final selection 


5 6.63±0.26 1.29 


5.34 


58.1 


x h > -1.8 


1 2.50±0.17 0.58 


1.92 


49.6 




Muon channel 158.0 pb~ 


i 




Preselection 


6441 6177 4871 


1239 


84.8 


tight lepton id. 


15 15.5 ±0.7 1.96 


13.6 


75.0 


final selection 


5 5.09 ±0.19 0.09 


5.00 


70.8 


x h > -1.74 


2 1.69±0.12 0.02 


1.67 


60.5 


Tau channel 158.0 pb 1 


Preselection 


7128 7091 4810 


2281 


95.8 


£ + £~qq 


21 20.4± 0.5 3.8 


16.6 


31.4 


final selection 


11 11.54±0.39 1.73 


9.81 


29.9 


C> 0.83 


0.77±0.03 0.03 


0.74 


18.1 


Missing energy channel 153.3 pb 1 


Anti 77 


14294 13623.2 10854.6 


2563.2 


84.3 


Preselection 


1183 1152.9 705.9 


430.9 


77.3 


C> 2.55 


27 27.8± 1.0 17.1 


10.1 


54.3 


C> AA 


4 6.0±0.22 3.1 


2.9 


33.9 




Four-jet channel 158.0 pb 


-i 




Preselection 


1730 1706.2 583 


1123 


87.1 


C> -1.0 


136 122.9± 1.1 26.8 


96.1 


63.3 


£> 0.28 


24 24.9± 0.2 7.2 


17.7 


45.9 


hA four-jet channel 158.0 pb 1 


Preselection 


1327 1274 318 


956 


85.9 


£> 2.0 


13 11.41± 0.34 5.19 


6.22 


65.5 


C> 3.05 


3 4.72± 0.15 2.11 


2.61 


55.0 



Table 1: Effect of the selection cuts on data, simulated background and simulated sig- 
nal events at y/s = 188.7 GeV. Efficiencies (in %) are given for the signal, ie. 
m H = 95 GeV/c 2 for the SM and m A = 80 GeV/c 2 , tan/3 = 20 for the MSSM. Within 
each channel, the last line gives the entries for the mass-plot, while the preceding line 
represent the inputs for the limit derivation. The quoted errors are statistical only. 
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m H Electron Muon Ht+t r+r Z Mis. Energy Four-jet 

(GeV/c 2 ) channel channel channel channel channel channel 

70.0 55.4^;? 68.3l{| 28.1 ± 3.0 32.0 ± 3.5 20.6 ± 2.0 52.2 ± 4.1 

75.0 56.81^ 71.2±i;| 28.3 ± 3.0 30.8 ± 3.3 32.3 ± 3.0 54.6 ± 4.2 

80.0 58.11^ 73.4±i;g 28.1 ± 3.0 31.5 ± 3.5 43.5 ± 4.0 58.7 ± 4.5 

85.0 57. Bt^ 72.2t\ A 4 27.6 ± 3.1 29.9 ± 3.2 52.0 ± 4.6 58.8 ± 4.5 

90.0 59.0^i 73.7t\i 26.6 ± 2.8 30.6 ± 3.2 57.1 ± 5.1 62.4 ± 4.7 

95.0 58.liy 70.81^ 25.9 ± 2.7 29.9 ± 3.1 54.3 ± 4.8 63.3 ± 4.7 

100.0 55.5^| 62.01^ 26.3 ± 2.8 28.2 ± 3.0 45.5 ± 4.1 56.7 ± 4.3 

Table 2: hZ channels: efficiencies (in %) of the selection at ^/s = 188.7 GeV as a function 
of the mass of the Higgs boson. The quoted errors include systematic uncertainties. 



tan 0= 20 


tan (3= 2 


m A 


Four-jet 


Tau 


m A 


Four-jet 


Tau 


(GeV/c 2 ) 


channel 


channel 


(GeV/c 2 ) 


channel 


channel 


70.0 


58.6 ±5.0 


32.4 ±3.4 


70.0 


51.8 ±4.4 


13.5 ±1.5 


75.0 


60.5 ±5.0 


33.1 ±3.4 


75.0 


54.7 ±4.5 


16.7 ± 1.8 


80.0 


65.5 ±5.3 


32.2 ±3.3 


80.0 


58.4 ±4.8 


25.8 ±2.7 


85.0 


64.7 ±5.3 


31.8 ±3.4 


85.0 


60.0 ±4.9 


33.7 ±3.6 


90.0 


60.6 ±5.0 




90.0 


61.1 ± 5.1 





Table 3: hA channels: efficiencies (in %) of the selection at y/s = 188.7 GeV as a 
function of the mass of the A boson for two values of tan (3(20 and 2). The quoted errors 
include systematic uncertainties. 
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the last line gives the result of a tighter selection. The efficiencies versus Higgs mass at 
the final selection level can be seen in Table |2] (SM channels) and Table |3] (MSSM chan- 
nels). The errors are obtained by summing the statistical and systematic uncertainties 
quadratically. 



8.2 The SM Higgs boson 



As an illustration, Figure [14] shows the distribution of the reconstructed Higgs boson 
mass found in the HZ channel after the tight selection (last lines in Table p]) for data, 
simulated background and signal events. The last cut in each channel has been chosen 
such that the signal-over-background ratio (for the reference mass) be almost equal for 
all channels (between .2 and .35). Within the 188.7 GeV data, the total number of events 
observed in all channels is 31, which is consistent with the 35.9 expected background. 
Since the two hypotheses (background only and background plus signal at 95 GeV/c 2 ) 
are almost indistinguishable, a possible signal at 90 GeV/c 2 has been superimposed in 
order to visualise our resolution in the mass of such a signal. 

We proceed to set a limit on the SM Higgs boson mass, combining these data with 
those taken at lower energies, namely 161,172 GeV | |14|| a nd 183 GeV The expected 
cross-sections and branching ratios are taken from ||2d|j2T|| , with the top mass set to 
175 GeV/c 2 . 

The confidence levels CL b , CL sb and CL S are computed as described in 0. CL b 
and CL sb are the confidence levels in these hypotheses (background only and signal + 
background), while CL S is conservatively taken as their ratio (CL sb /CL b ). 

In the presence of a sizeable Higgs signal the value of the observed CL b (top of Fig- 
ure |ll|) would approach one, because it measures the fraction of background-only ex- 
periments which are less signal-like than the observation. On the contrary here, the 
observation agrees well with the expectation (background only). Furthermore the curve 
for the signal hypothesis shows that the expected 5a discovery limit (horizontal line at 
1 — CL b = 5.7 x 10~ 7 ) is at 88.6 GeV/c 2 . The confidence level in the signal is shown in 
Fig. [15] (bottom). The observed 95% CL lower limit on the mass is mn > 94.6 GeV/c 2 , 
while the expected mean is 94.4 GeV/c 2 and the expected median (50% exclusion poten- 
tial) is 95.3 GeV/c 2 . If errors had not been allowed for, the observed (expected) limit 
would have been increased by 0.2 GeV/c 2 (0.4 GeV/c 2 ). 

The effective A% 2 (—2 A In C) with which the SM Higgs is excluded is shown in Fig. [16. 
In the event of a discovery the A% 2 would be negative, and could be used to extract the 



mass and its error, as can be seen on the bottom plot of Fig. |T6" 



Finally the data can be used to set 95% CL upper bounds on the HZZ coupling in 
non-standard models which assume that the Higgs boson decay properties are identical 
to those in the SM but the cross-section may be different. Figure |17] shows the excluded 
cross-section function of the test mass. 



8.3 Neutral Higgs bosons in the MSSM 

The results in the hZ and hA processes are combined with the same statistical method 
as for the SM, using also earlier results at LEP2 energies [|I],|T^.^]. 

In the MSSM, at tree level, the production cross-sections and the Higgs branching 
fractions depend on two free parameters, tan/3 and one Higgs boson mass, or, alterna- 
tively, two Higgs boson masses, eg and rrv The properties of the MSSM Higgs 
bosons are modified by radiative corrections which introduce additional parameters: the 
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mass of the top quark, the Higgs mixing parameter, /i, the common sfermion mass term 
at the EW scale, Ms, the common SU(2) gaugino mass term [] at the EW scale, M 2 , 
and the common squark trilinear coupling at the EW scale, A. The interpretation of the 
experimental results depends on the values assumed for these parameters as well as on 
the order of the calculated radiative corrections. 

The results described hereafter rely on leading-order two-loop calculations of the ra- 
diative corrections in the renormalization group approach , with recent modifications 



(about top threshold and gluino two-loop corrections) that make the computations agree 
with fully diagrammatic two-loop calculations |24| . After these improvements, the bench- 
mark prescriptions for the parameters beyond tree-level have also been refined, leading 
to two extreme scenarii for the theoretical upper bound on m n as a function of tan f3 p5| 



which differ only by the value of X t = A — /xcot/3, the parameter which controls the 
mixing in the stop sector. In the following, we adopt these new prescriptions which cor- 
respond to: 175 GeV/c 2 for the top mass, 1 TeV/c 2 for M s , 200 GeV/c 2 for M 2 and 
-200 GeV/c 2 for \i. Two values have been considered for the mixing in the stop sector: 
X t = ^/QMs, which defines the so-called m™ ax scenario, and X t = 0, which defines the 
no mixing scenario. Then a scan is made over the MSSM parameters tan/? and , in 
the wa range [| of 20 GeV/c 2 - 1 TeV/c 2 , and tan/3 between 0.5 and 50. At each point 
of the parameter space, the hZ and hA cross-sections and the Higgs branching fractions 



are computed with the HZHA03 |T(J program. 

The signal expectations in each channel are derived from the cross-sections, the ex- 
perimental luminosity and the efficiencies. A correction is applied to account for differing 
branching fractions of the Higgs bosons into bb and r + r~ between the input point and 
the simulation (e.g. for the hZ process, the simulation is done in the SM framework). 
For the hA channels, as there can be a difference between the masses of the h and A 
bosons at low tan /?, the set of hA efficiencies obtained from the simulation at tan (3 = 20 
is applied at all points with tan/3 above 2.5, while the set of efficiencies derived from the 
tan f3 — 2 simulation is applied below. The same holds for the discriminant information. 
The signal expectations, expected backgrounds and numbers of candidates enter in the 
computation of the observed confidence level in the signal hypothesis at the input point, 
CL S . The expected confidence level in the signal hypothesis is also derived at each point. 
As there is a large overlap in the background selected by the two four-jet channels, only 
one channel is selected at each input point, on the basis of the best signal over back- 
ground ratio. This ensures that the channels which are combined in the confidence level 
computations are independent. 

The results translate into regions of the MSSM parameter space excluded at 95% CL 



or more. The excluded regions are presented in the (m n , tan j3) plane in Fig. [18], in the 



toa , tan/?) plane in Fig. |1^ and in the (m\ , m^) plane in Fig. |20[. As illustrated in the 
latter, there is a small region of the parameter space where the decay h-^AA opens, in 
which case it supplants the h— >bb decay. But, due to the high luminosity collected at 
188.7 GeV, the results in the h— >bb channel alone cover most of the area which remained 
unexcluded at 183 GeV in this region [I]]. 

Finally, the results shown in Figs. ITSL [T9] and |20 establish 95% CL lower limits on 



and rriA , whatever the assumption on the mixing in the stop sector and for all values of 
tan/? greater than or equal to 0.6: 

m h > 82.6 GeV/c 2 m A > 84.1 GeV/c 2 . 



1 The U(l) and SU(3) gaugino mass terms at the EW scale, Mi and M3, are assumed to be related to M2 through the 
GUT relations Mi = (5/3)tan 2 w M 2 and M 3 = (a s /a)sin 2 w M 2 . 

2 The region m\ below 20 GeV/c 2 would need LEP1 results which are not yet available in the format required by the 
statistical procedure. 
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The expected limits are 81.3 GeV/c 2 in m h and 82.3 GeV/c 2 in . In the low tan/3 
region, in the no mixing case, all values of up to 1 TeV/c 2 are excluded, providing an 
excluded range in tan/3 between 0.6 and 2.2, in agreement with the expected excluded 
range. On the other hand, no limit can be set on tan/? in the m™ ax scenario (see Fig. [3~Q] ) . 



8.4 Interpretation in a general Two Higgs Doublet Model 

These results can also be translated to the framework of a general Two Higgs Doublet 
Model (2HDM) with one assumption (the decay of both h and A is dominated by bb 
and/or cc final states) and two options: CP conserving or CP violating. 

In the CP-conserving two Higgs doublet model, the h and H bosons are mixtures of the 
real parts of the neutral Higgs fields, while A derives from the imaginary components not 
absorbed by the Z. The coupling strengths are: C^z = sin ( a ~ 0) an d C^a = cos (a— (3). 
These couplings clearly indicate the complementarity of the two processes. Besides, if 
one of them is experimentally out of reach for a given set of masses, no exclusion is 
possible since mixing angles could always be such that the other process is suppressed 
below detectability. 

The exclusion plot (left of Figure |2T[) is obtained in the following way: for each pair of 
mh and values the number of expected events for the channels hA and hZ is calculated 
using the cross-sections, integrated luminosities, branching ratios and efficiencies quoted 
in this paper. This number depends obviously on the factors C^z and C^A • Thus a 
minimisation of the confidence level with respect to these factors has been performed, 
taking into account the sum rule C?^ + ^\p. = ^' 

It should be noted that in 2HDM the branching ratios of A and h into bb are propor- 
tional to tan (3 and sin a / cos (3 respectively. Imposing the condition | sin a \ > cos /3, which 
is barely restrictive for medium or large values of tan /3, leads to a dominant coupling to 
b quarks for both h and A bosons, (Zone I in Figure pip . 

The case of non-b decays has also been studied, using the selections of this paper 
except those referring to b-tagging. This takes care of a scenario with tan/3< 1 which 



would allow a dominant decay of the Higgs boson into cc |26|. This region (named Zone 



II in Figure |2l| ) occurs for values of the a and (3 angles such that a ~ (3 ~ 0, resulting in 
a very restrictive and particular parameter set of the 2HDM which represents the most 
pessimistic scenario in this kind of search. Any other situation will lead to an intermediate 
excluded region, as for example when sin a = but tan (3 > 1 which implies qq bb decays 
of the hA signal. 

CP violation in the SUSY sector is an open possibility and may even be necessary 
in the electroweak baryogenesis scenario (see for instance p7 and references therein). 



The violation leads to three neutral Higgs bosons (noted hi, I12 and li3, sorted by mass) 
with undefined CP properties. It is shown in p8| that the previous sum rule, valid in 
the CP-conserving case, can be extended to the CP-violating model, giving the relation: 
ChiZ + Ch 2 z + Chih 2 = 1' which together with the results described in this paper allows the 
minimum number of events for each point (m^, mh 2 ) to be calculated, with an extended 
procedure similar to the one used previously, leading to the excluded region at 95% CL 
shown in the lower plot of Figure |2T] (only the conditions for Zone I were used in this 
case) . 

This analysis shows that it is possible to exclude a large region of Higgs boson masses 
even when relaxing the standard assumptions made in the MSSM scenario. 
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9 Conclusions 

The 158 pb _1 of data taken by DELPHI at 188.7 GeV, combined with our lower energy 
data, sets the lower limit at 95% CL on the mass of the Standard Model Higgs boson at: 

m H > 94.6 GeV/c 2 . 

The MSSM studies described above give for all values of tan (3 above 0.6, and assuming 
m A >20GeV/c 2 : 

m h > 82.6 GeV/c 2 m A > 84.1 GeV/c 2 . 

Beyond the results described above, DELPHI performed a more complete scan of the 
parameters of the MSSM: this analysis, described in the addendum, shows clearly the 
robustness of the limits obtained in the benchmark scenarios. 

Other LEP experiments, using their data sets collected concurrently with the ones 



used in this work, have reported similar results [29,30] 
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Figure 1: Distributions of the combined b-tagging variable x\>, in data (points) and sim- 
ulation (histogram) . The contribution of udsc-quarks is shown as the dark histogram. 
Bottom: the ratio of the tagging rates in the data and the simulation as a function of the 
cut in the b-tagging variable. 
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Figure 2: He + e~channel: distributions of some analysis variables as described in the 
text, at the preselection level. The plots on the left-hand side show a comparison between 
188.7 GeV data (dots) and simulated background events (solid line) normalised to the 
experimental luminosity. The dark grey areas represent the contribution of the qq(7) 
background and the light grey area the e + e~qq contribution. The expected normalised 
distributions of the same variables for a signal at 95 GeV/c 2 are represented on the 
right-hand side. Note the different y-scales. 
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Figure 3: He + e~channel: distributions of the global b-tagging variable versus the fitted 
recoil mass for data, background and simulated signal events with m H = 95 GeV/c 2 at 
188.7 GeV. These distributions are the inputs for the Confidence Levels computation. 
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Figure 4: H/i + yU~ channel: distributions of some analysis variables as described in the 
text, at the preselection level (a to f) or after the lepton pair selection (g and h). The 
dark grey area represents the contribution of the qq(7) background and the light grey area 
the 4 fermions contribution. The expected normalised distribution for the signal (with 
mn =95 GeV/c 2 ) is represented on the right-hand column. 
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Figure 6: Yivv channel: distributions of the main analysis variables as described in the 
text, at the preselection level. The points with error bars represent the data. The left 
hand side histograms show the different backgrounds and the right hand side histograms 
show the signal distributions for mn = 95 GeV/c 2 . 
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Figure 7: Hz/z/ channel: Distributions of the discriminant variable for the ex- 
pected background events (full histograms), data (points) and a Higgs signal with mn = 
95 GeV/c 2 (dashed histogram) are shown. The normalisation for the Higgs signal is 100 
times the expectation. 




Figure 8: Hz/z/ channel: Efficiency versus total background as a result of varying the cut 
on the Likelihood. The contribution of different backgrounds are shown separately. Points 
with error bars are the data, and the thick black line shows the total expected background. 
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Figure 9: Hqq channel: distributions of the main analysis variables at the preselection 
level. The points with error bars represent the data. The left hand side histograms show 
the different backgrounds and the right hand side ones show the expected signal distribu- 
tions for m-H =95 GeV/c 2 . 
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Figure 10: Hqq channel: The reconstructed mass at the preselection level. Note that the 
signal corresponds to 400 times the luminosity for clarity. 
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Figure 11: Hqq channel: Efficiency versus total background as a result of varying the cut 
on the Likelihood. The contribution of different backgrounds are shown separately. Points 
with error bars are the data and the upmost black line the total expected background. 
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Figure 12: hA channel: Efficiency versus total background as a result of varying the cut 
on the Likelihood. The contribution of different backgrounds are shown separately. Points 
with error bars are the data and the upmost black line the total expected background. 
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Figure 13: hA hadronic channel: Mass distribution for the sum of the di-jets in data and 
MC at the end of the hA J^-jet analysis at the 55% efficiency level. The signal, with = 
80 GeV/c 2 and tan/3 = 20 is normalised to the data. 
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Figure 14: Final distribution of the reconstructed Higgs boson mass when combining all 
hZ analyses at 182.7 and 188. 7 GeV. Data are compared with background expectations. 
The expected spectrum, with the correct rate, from a signal at 90 GeV/c 2 is also shown 
added to the background contributions, as the dotted histogram. 
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Figure 15: Confidence levels as a function of the SM Higgs boson mass. Curves are shown 
for the expected (dashed) and observed (solid) confidences and the bands correspond to the 
68.3% and 95% confidence intervals. Top: Confidence level in the background hypothesis. 
Bottom: Confidence level in the signal hypothesis. The intersections of the horizontal line 
at 5% with the curves define the expected and observed 95% CL lower limits on the Higgs 
boson mass. 
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Figure 16: The effective Ax 2 with which each SM Higgs mass is excluded (solid) and the 
expected value of the same (dashed). Top: Expectation in case of background only; the 
dark/light bands correspond to the 68.3% / 95% confidence intervals. Bottom: the bands 
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Figure 17: 95% CL excluded cross-sections as a function of the Higgs boson mass compared 
with the SM expectation. 
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Figure 18: Regions in the (m^, tan/3) plane excluded at 95% CL by the searches in 
the hZ and hA channels up to \/s = 189 GeV. Two extreme hypotheses for the mix- 
ing in the stop sector are considered. The regions not allowed by the MSSM model for 
mtop = 175 GeV/c 2 , M SUS y = 1 TeV/c 2 ; M 2 = -// = 200 GeV/c 2 andm A < 1 MeV/c 2 or 
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Figure 19: Regions in the (m^, tan (3) plane excluded at 95% CL by the searches in the 
hZ and hA channels up to = 189 GeV. Two extreme hypotheses for the mixing in the 
stop sector are presented. 



32 



M « 250 
> 

Qi 

S 225 

E 

200 
175 


F>FT PHT 

1 / 1 V 1 ^ I 111 

Vs from 130 to 189 GeV 

; m top = 175GeV/c 2 
: M susy = lTeV/c 2 

M 2 = -n = 200 GeV/c 

tan (3 > 0.5 
: No mixing 1 






150 


— observed limit / 
expected limit / 






125 








100 


1 theoretically forbidden / ,y 


75 




50 


/ excluded / 
1 / \ theoretically forbidden 


25 


^ , , i S//Y/A , , a , , 



20 40 60 80 100 120 140 



m h (GeV/c 2 ) 











« 250 


DELPHI 






> 

£ 225 


Vs from 130 to 189 GeV 






< 

s 

200 


; m top = 175GeV/c 2 
i M susy = lTeV/c 2 

M 2 = -\i = 200 GeV/c 2 / 










175 


tan (3 > 0.5 / 

: m.max scenario / 

/ 






150 


— observed limit / 
expected limit / 






125 








100 








75 


theoretically forbidden / / 
/ excluded / 






50 


1 / theoretically forbidden 


25 










20 40 60 80 100 120 140 



m h (GeV/c 2 ) 

Figure 20: Regions in the (m^ , m^) plane excluded at 95% CL by the searches in the 
hZ and hA channels up to y/s = 189 GeV (in light grey). Two extreme hypotheses for the 
mixing in the stop sector are presented. The regions not allowed by the MSSM model for 
m top = 175 GeV/c 2 , M SUSY = 1 TeV/c 2 , M 2 = -fi = 200 GeV/c 2 and m A < 1 MeV/c 2 or 
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Figure 21: Excluded regions at the 95% CL a) in the CP-conserving model (with dominant 
b-decays /zone I) and dominant non-b decays (zone II)) and b) in the CP-violating 



